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Abstract. With the new generation of instruments for 
Cosmic Microwave Background (CMB) observations aim- 
ing at an accuracy level of a few percent in the measure- 
ment of the angular power spectrum of the anisotropics, 
the study of the contributions due to secondary effects 
has gained impetus. Furthermore, a reinvestigation of the 
main secondary effects is crucial in order to predict and 
quantify their effects on the CMB and the errors that they 
induce in the measurements. 

In this paper, we investigate the contribution, to the 
CMB, of secondary anisotropics induced by the transverse 
motions of clusters of galaxies. This effect is similar to the 
Kaiser-Stebbins effect. In order to address this problem, 
we model the gravitational potential well of an individ- 
ual structure using the Navarro, Frenk & White profile. 
We generalise the effect of one structure to a population 
of objects predicted using the Press-Schechter formalism. 
We simulate maps of these secondary fluctuations, com- 
pute the angular power spectrum and derive the average 
contributions for three cosmological models. We then in- 
vestigate a simple method to separate this new contribu- 
tion from the primary anisotropics and from the main sec- 
ondary effect, the Sunyaev-Zel'dovich kinetic effect from 
the lensing clusters. 

Key words: Cosmology: cosmic microwave background - 
gravitational lensing - secondary fluctuations - clusters of 
galaxies 



1. Introduction 

During the next decade, several experiments are planned 
to observe the Cosmic Microwave Background (CMB) and 
measure its temperature fluctuations (Planck surveyor, 
Map, Boomerang, ...). Their challenge is to measure the 
small scales anisotropics of the CMB (a few arcminutes 
up to ten degrees scale) with sensitivities better by a 
factor 10 than the COBE satellite (Smoot et al. 1992). 
These high sensitivity and resolution measurements will 
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tightly constrain the value of the main cosmological pa- 
rameters (Kamionkowski et al. 1994). However, the con- 
straints can only be set if we are able to effectively measure 
the primary temperature fluctuations. These fluctuations, 
present at recombination, give an insight into the early 
universe since they are directly related to the initial den- 
sity perturbations which are the progenitors to the cosmic 
structures (galaxies and galaxies clusters) in the present 
universe; but which are first and foremost the relics of the 
very early initial conditions of the universe. 
Between recombination and the present time, the CMB 
photons could have undergone various interactions with 
the matter and structures present along their lines of sight. 
Some of these interactions can induce additional tem- 
perature fluctuations called, secondary anisotropies be- 
cause they are generated after the recombination. Along 
a line of sight, one measures temperature fluctuations 
which are the superposition of the primary and secondary 
anisotropies. As a result, and in the context of the future 
CMB experiments, accurate analysis of the data will be 
needed in order to account for the foreground contribu- 
tions due to the secondary fluctuations. Photon-matter 
interactions between recombination and the present time 
are due to the presence of ionised matter or to variations 
of the gravitational potential wells along the lines of sight. 

The CMB photons interact with the ionised matter 
mainly through Compton interactions. In fact, after re- 
combination the universe could have been re-ionised glob- 
ally or locally. Global early re-ionisation has been widely 
studied (see Dodelson & Jubas 1995 for a recent re- 
view and references therein). Its main effect is to either 
smooth or wipe out some of the primary anisotropies; 
but the interactions of the photons with the matter in 
a fully ionised universe can also give rise to secondary 
anisotropies through the Vishniac effect (Vishniac 1987). 
This second order effect has maximum amplitudes for 
a very early re-ionisation. The case of a late inhomoge- 
neous re-ionisation and its imprints on the CMB fluctu- 
ations has been investigated (Aghanim et al. 1996) and 
found to be rather important. In this case, the secondary 
anisotropies are due to the bulk motion of ionised clouds 
with respect to the CMB frame. When the re-ionisation 
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is localised in hot ionised intra-cluster media the pho- 
tons interact with the free electrons. The inverse Compton 
scattering between photons and electrons leads to the so- 
called Sunyaev-Zel'dovich (hereafter SZ) effect (Sunyaev 
& Zel'dovichl972, 1980). The Compton distortion due to 
the motion of the electrons in the gas is called the ther- 
mal SZ effect. The kinetic SZ effect is a Doppler distor- 
tion due to the peculiar bulk motion of the cluster with 
respect to the Hubble flow. The SZ thermal effect has the 
unique property of depressing the CMB brightness in the 
Raylcigh- Jeans region and increasing its brightness above 
a frequency of about 219 GHz. This frequency dependence 
makes it rather easy to observe and separate from the ki- 
netic SZ effect. In fact, the latter has a black body spec- 
trum which makes the spectral confusion between kinetic 
SZ and primary fluctuations a serious problem. The SZ 
effect has been widely studied for individual clusters and 
for populations of clusters. For full reviews on the subject 
we refer the reader to two major articles: Rephaeli 1995 
and Birkinshaw 1997. These investigations have clearly 
shown that the SZ effect in clusters of galaxies provides a 
powerful tool for cosmology through measurements of the 
Hubble constant, the radial peculiar velocity of clusters 
and consequently the large scale velocity fields. 

Besides the interactions with the ionised matter, some 
secondary effects arise when the CMB photons traverse a 
varying gravitational potential well. In fact, if the grav- 
itational potential well crossed by the photons evolves 
between the time they enter the well and the time they 
leave it, the delay between entrance and exit is equiva- 
lent to a shift in frequency, which induces a temperature 
anisotropy on the CMB. This effect was first studied by 
Rees & Sciama (1968) for a potential well growing under 
its own gravity. Numerous authors have investigated the 
potential variations due to collapsing objects and their ef- 
fect on the CMB (Kaiser 1982, Nottale 1984, Martinez- 
Gonzalez, Sanz & Silk 1990, Seljak 1996). Similarly, a 
gravitational potential well moving across the line of sight 
is equivalent to a varying potential and will thus imprint 
secondary fluctuations on the CMB. This effect was first 
studied for one cluster of galaxies by Birkinshaw & Gull 
(1983) (Sect. 2). Kaiser & Stebbins (1984) and Bouchet, 
Bennett & Stebbins (1988) investigated a similar effect 
for moving cosmic strings. Recent work (Tuluic & Laguna 
1995, Tuluie, Laguna & Anninos 1996) based on N-body 
simulations has pointed out this effect in a study of the 
effect of varying potential on rather large angular scales 
(~ 1°). A discussion of some of these results and a com- 
parison with ours will follow in the next sections. 

In this paper, following the formalism of Birkinshaw 
& Gull (1983) and Birkinshaw (1989), we investigate the 
contribution of secondary anisotropies due to a popula- 
tion of collapsed objects moving across the line of sight, 
these objects range from small groups to rich clusters in 
scale (10 13 to 10 15 Mq). In section 2., we first study in 
detail the case of a unique collapsed structure. We use a 



structure model to compute in particular the deflection 
angle and derive the spatial signature of the moving lens 
effect. We then account (Section 3.) for the contribution, 
to the primordial cosmological signal, of the whole popu- 
lation of collapsed objects using predicted counts and we 
simulate maps of these secondary anisotropies. In section 
4., we analyse the simulated maps and present our results. 
We give our conclusions in section 5. 

2. Formalism for an individual moving structure 

One of the first studies of the photon-gravitational po- 
tential well interactions is related to the Sachs- Wolfe ef- 
fect (Sachs & Wolfe 1967). At the recombination time 
(z ~ 1100) the photons and matter decouple while they 
are in potential wells; the photons are redshifted when 
they leave the potential wells. This generates the large 
angular scale temperature fluctuations. 
Other authors have investigated the effect of time varying 
potentials on the CMB photons after the recombination, 
namely the Rees-Sciama effect (Rees & Sciama 1968). If 
the potential well crossed by the photons evolves between 
the time they enter and their exit, the extra-time delay 
they suffer changes the temperature of the CMB and in- 
duces an additional anisotropy. The variation of the po- 
tential well can have an "intrinsic" or a "kinetic" origin. 
The first case describes the evolution with respect to the 
background density distribution. The second case is re- 
lated to the bulk motion of a gravitational potential well 
across the line of sight which mimics a time variation of 
the potential. Photons crossing the leading edge of a struc- 
ture will be redshifted because of the increasing depth of 
the potential well during their crossing time; while pho- 
tons crossing the trailing edge of the same structure are 
blueshifted. This results in a characteristic spatial signa- 
ture for the induced anisotropy: a hot-cold temperature 
spot. 

The specific effect of a moving cluster across the sky was 
first studied by Birkinshaw & Gull (1983) (correction to 
this paper was made in Birkinshaw 1989) and it was in- 
voked as a method to measure the transverse velocity of 
massive clusters of galaxies. These authors found that the 
transverse motion of a cluster across the line of sight in- 
duces a frequency shift given by: 

Ais 

= /?7 sin a cos (p5(b). (1) 

Here, (i is the peculiar velocity in units of the speed of light 
(f3 = v/c), 7 is the Lorenz factor (7 = (1 — /3 2 ) -1 / 2 ), a and 
4> are respectively the angle between the peculiar velocity 
v and the line of sight of the observer and the azimuthal 
angle in the plane of the sky, and (5(b) is the deflection 
angle due to the gravitational lensing by the cluster at a 
distance equal to the impact parameter b. This frequency 
shift induces a brightness variation which in turn can be 
expressed as a secondary temperature fluctuation ST/T. 
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Fig. 1. Characteristic spatial signature of a temperature fluc- 
tuation due to a moving lens with mass M = 10 15 Mq and 
velocity v = 600 km/s. 

In their paper, Birkinshaw & Gull derived an expression 
for 8T/T in the Rayleigh- Jeans regime, with some spe- 
cific assumptions on the gravitational potential well asso- 
ciated with the cluster. They assumed that the matter in 
the galaxy cluster was homogeneously distributed in an 
isothermal sphere of radius i?, where R is the characteris- 
tic scale of the cluster. 

In our paper, we basically follow the same formalism as 
Birkinshaw & Gull's using the corrected expression from 
Birkinshaw 1989. We compute the gravitational deflection 
angle at the impact parameter 5(h), the corresponding fre- 
quency shift and then derive the associated temperature 
fluctuation. The main difference between our approach in 
this section and the previous work concerns the physical 
hypothesis that we adopt to describe the distribution of 
matter in the structures. In fact, in order to derive the de- 
flection angle, we find the homogeneous isothermal distri- 
bution a too simple and rather unrealistic hypothesis and 
choose another more realistic description. For the struc- 
tures such as those we are interested in (clusters down 
to small groups), almost all the mass is "made" of dark 
matter. In order to study the gravitational lensing of a 
structure properly, one has to model the gravitational po- 
tential well using the best possible knowledge for the dark 
matter distribution. The corrections, due to the more ac- 
curate profile distribution that we introduce, will not alter 
the maximum amplitude of an individual moving lens ef- 
fect since it is associated with the central part of the lens. 
However, when dealing with some average signal coming 
from these secondary anisotropies, the contribution from 
the outskirts of the structures appears important and thus 
a detailed model of the matter profile is needed. 
In view of the numerous recent studies on the formation 
of dark matter halos, which are the formation sites for 
the individual structures such as clusters of galaxies, we 
now have a rather precise idea of their formation and den- 



sity profiles. Specifically, the results of Navarro, Frenk & 
White (1996, 1997) are particularly important. In fact, 
these authors have used N-body simulations to investigate 
the structure of dark matter halos in hierarchical cosmogo- 
nies; their results put stringent constraints on the dark 
matter profiles. Over about four orders of magnitudes in 
mass (ranging from the masses of dwarf galaxy halos to 
those of rich clusters of galaxies) , they found that the den- 
sity profiles can be fitted over two decades in radius by a 
"universal" law (hereafter NFW profile) which seems to 
be the best description of the structure of dark matter 
halos (Huss, Jain & Steinmetz 1997). The NFW profile is 
given by: 



p{r) 



Pcrit $c 



(r/r s ){l + r/r s ) 2 



(2) 



where r s = 7"2oo/c is the scale radius of the halo, S c its 
characteristic overdensity, p cr it is the critical density of 
the universe and c is a dimensionless parameter called 
the concentration. The radius r2oo is the radius of the 
sphere where the mean density is 200 x p cr n. This is 
what we refer to as a virialised object of mass M200 = 
200p CI , t (47r/3)r| 00 . 

In addition to the fact that the shape is independent of 
the halo mass over a wide range, the NFW profile is also 
independent of the cosmological model. The cosmological 
model intervenes essentially in the formation epoch of the 
dark matter halo and therefore in the parameters of the 
profile, namely c, r s and S c . 

Using the density profile, one can compute the deflec- 
tion angle at the impact parameter which gives the shape 
of the pattern and the amplitude of the induced secondary 
anisotropy. In our work, we compute the deflection angle 
following the formalism of Blandford & Kochanek (1987), 
which is given by the expression: 



2 ZT Vr / $ ( r '^' 



(3) 



here, the integral is performed over the length element dl 
along the line of sight. Di s and D os are respectively the 
distances between lens and source and the observer and 
source. In the redshift range of the considered structures 
(z < 1.5), the distance ratios Di s /D os range between 1 
and 0.68 for the standard CDM model, between 1 and 
0.53 for the open CDM and between 1 and 0.74 for the 
lambda CDM model. These cosmological models will be 
defined in the next section. In Eq. [| r is the position of 
the structure and $(r, I) is the associated gravitational 
potential. In order to get an analytic expression of the 
deflection angle and hence of the anisotropy, we used a 
density profile which gives a good approximation to the 
NFW density profile (Eq. ||), in the central part of the 
structure. This density profile is given by: 



p(r) = Pent S c 



exp 



(4) 
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The fitted profile leads to a diverging mass at large radii 
and we therefore introduce a cut-off radius R max to the 
integral. This cut-off should correspond to some physical 
size of the structure. With regard to the different values of 
the concentration c, we set R m ax = 8r s which is in most 
cases equivalent to Rmax — ?"200,*- e -; close to the virial ra- 
dius. The integral giving the deflection angle is performed 
on the interval [—Rmax, Rmax]- For R max = 8r s , our fit 
gives a mass which is about 20% lower than the mass de- 
rived from NFW profile. This difference is larger for larger 
Rmax, and for R m ax = 10r s we find that the mass is about 
33% lower. However, the larger radii the temperature fluc- 
tuations are at the 10 -8 level. On the other hand, the 
Hernquist (1990) profile is also in agreement with the re- 
sults of N-body simulations. Indeed, both NFW and Hern- 
quist profiles have a similar dependence in the central part 
of the structure but differ at large radii where the NFW 
profile is proportional to r~ 3 and the Hernquist profile 
varies as r -4 . However, the amplitude of the anisotropy 
at large radii is very small and the results that we obtain 
does are not sensitive to the cut-off. 

Given the peculiar velocity of the structure and its 
density profile, we can calculate the deflection angle (Eq. 
||). Then one can determine the relative variation in fre- 
quency, 8v/v 7 using equation [j] and thus evaluate the sec- 
ondary distortion induced by a specific structure moving 
across the sky. We find that individual massive structures 
(rich galaxy clusters) produce anisotropies ranging be- 
tween a few 10 -6 to 10 -5 ; but within a wider range of 
masses the amplitudes are smaller and these values are 
only upper limits for the moving lens effect. 



3. Generalisation to a sample of structures 

Future CMB (space and balloon born) experiments will 
measure the temperature fluctuations with very high ac- 
curacy (10~ 6 ) at small angular scales. In our attempts to 
foresee what the CMB maps would look like and what 
would be the spurious contributions due to the various as- 
trophysical foregrounds, we investigate the generalisation 
of the computations made above to a sample of structures. 
This is done in order to address the questions of the cu- 
mulative effect and contamination to the CMB. 

Some work has already been done by Tuluie & Laguna 
1995 and Tuluie, Laguna & Anninos 1996 who pointed 
out the moving lens effect in their study of the vary- 
ing potential effects on the CMB. In their study, they 
used N-body simulations to evolve the matter inhomo- 
geneities, from the decoupling time until the present, in 
which they propagated CMB photons. They have esti- 
mated the anisotropies generated by three sources of time- 
variations of the potential: intrinsic changes in the gravi- 
tational potential, decaying potential effect from the evo- 
lution of gravitational potential in JJo ^ 1 models, and pe- 
culiar bulk motions of the structures across the sky. They 
evaluated the contribution of the latter effect for rather 



large angular scales (~ 1°) due to the lack of numerical 
resolution (about 2/i _1 Mpc) and gave estimates of the 
power spectrum of these effects. 

With another approach, we make a similar analysis 
in the case of the moving lens effect extended to angular 
scales down to a few tens of arcseconds. We also simulate 
attempts at the detection and subtraction of the moving 
lens effect. Our approach is quite different from that of Tu- 
luie, Laguna & Anninos, in that it is semi-empirical and 
apply the formalism developed for an individual structure 
(Sect. 2.) to each object from a sample of structures. The 
predicted number of objects in the sample being derived 
from the Press-Schechter formalism for the structure for- 
mation (Press & Schechter 1974). 

3.1. Predicted population of collapsed objects 

An estimate of the cumulative effect of the moving lenses 
requires a knowledge of the number of objects of a given 
mass that will contribute to the total effect at a given 
epoch. We assume that this number is accurately pre- 
dicted by the abundance of collapsed dark matter halos 
as a function of their masses and redshifts, as derived 
using the Press-Schechter formalism. This approach was 
used in a previous paper (Aghanim et al. 1997) which pre- 
dicted the SZ contribution to the CMB signal in a stan- 
dard CDM model. In addition to the "traditional" stan- 
dard Cold Dark Matter (CDM) model (fi = 1), in this 
paper we also address the question of a generalised mov- 
ing lens effect in other cosmological models. We extend 
the Press-Schechter formalism to an open CDM model 
(OCDM) with no cosmological constant (Slo = 0.3), and 
also a flat universe with a non zero cosmological constant 
(ACDM model) (fi = 0.3 and A = 0.7). Here VL is the 
density parameter, A is the cosmological constant given 
in units of 3Hq and Ho is the Hubble constant. We take 
Hq = lOO/i km/s/Mpc, and assume h = 0.5 throughout 
the paper. 

In any case, the general analytic expression for the num- 
ber density of spherical collapsed halos in the mass range 
[M, M + dM] can be written as (Lacey & Cole 1993): 



dn(M, z) 
dM 



2 p(z) din cr(M) S c0 (z) 
7r M 5 " din M o-(M) 



exp 



2a 2 (Af) 



(5) 



where ~p~(z) is the mean background density at redshift z 
and S c o(z) is the overdensity of a linearly evolving struc- 
ture. The mass variance a 2 (M) of the fluctuation spec- 
trum, filtered on mass scale Af, is related to the linear 
power spectrum of the initial density fluctuations P{k) 
through: 



ct 2 (M) 



1 

2^2 



k 2 P{k)W 2 {kR) dk, 
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where W is the Fourier transform of the window function 
over which the variance is smoothed (Peebles 1980) and R 
is the scale associated with mass M. In the assumption of 
a scale-free initial power spectrum with spectral index n, 
the variance on mass scale M can be expressed in terms of 
as, the rms density fluctuation in sphere of 8/i _1 Mpc size. 
The relationship between these two quantities is given by 
(Mathiesen & Evrard 1997): 

a(M) = (1.19no) a CT8M _Qt , 

with a = (n + 3) /6. It has been shown that as varies with 
the cosmological model and in particular with the density 
parameter f^o- A general empirical fitting function (og = 
AQq b ) was derived from a power spectrum normalisation 
to the cluster abundance with a rather good agreement in 
the values of the parameters A and B (White, Efstathiou 
& Frenk 1993, Eke, Cole & Frenk 1996, Viana & Liddle 
1996). In our work, we use the "best fitting values" from 
Viana & Liddle (1996) which are A = 0.6 and B = 0.36 + 
O.310o - 0.28^0 for an open CDM universe (fl < 1 and 
A = 0) or B = 0.59 - 0.16fi + 0.06Slg for a flat universe 
with a non zero cosmological constant (f^o + A = 1). We 
use n = — 1 for the spectral index in the cluster mass 
regime which is the theoretically predicted value. Some 
local constraints on the temperature abundance of clusters 
favour n = —2 (Henry & Arnaud 1991, Oukbir, Bartlett 
& Blanchard 1997) but we did not investigate this case. 

3.2. Peculiar velocities 

On the scale of clusters of galaxies, typically 8/i _1 Mpc, 
one can assume that the density fluctuations are in the 
linear regime. Therefore the fluctuations are closely re- 
lated to the initial conditions from which the structures 
arise. In fact, in the assumption of an isotropic Gaussian 
distribution of the initial density perturbations, the initial 
power spectrum P(k) gives a complete description of the 
velocity field through the three-dimensional rms velocity 
{vrms) predicted by the linear gravitational instability for 
an irrotational field at a given scale R (Peebles 1993). This 
velocity is given by: 



a(t)Hf(n,A) 



/ P(k)W 2 (kR)dk 



-1 1/2 



(6) 



where a(t) is the expansion parameter, the Hubble con- 
stant H and the density parameter Q, vary with time 
(Caroll, Press & Turner 1992). The function /(fi, A) is ac- 
curately approximated by /(fi, A) = fl os (Peebles 1980) 
even if there is a non zero cosmological constant (Lahav 
et al. 1991). Furthermore, under the assumptions of lin- 
ear regime and Gaussian distribution of the density fluc- 
tuations, the structures move with respect to the global 
Hubble flow with peculiar velocities following a Gaussian 
distribution f(y) — 

described by v r 



Tom exp( -2^r) which is full y 

This prediction is in agreement with 



numerical simulations (Bahcall et al. 1994, Moscardini et 
al. 1996). 

The present observational status of peculiar cluster 
velocities puts few constraints on the cosmological mod- 
els. Results from the Hudson (1994) sample using D„-cr 
and IRTF distance estimators give respectively v rms = 
688 ± 82 and 646 ± 120 km/s, a composite sample gives 
Vrms — 725±60 km/s (Moscardini et al 1996). Giovanelli's 

(1996) sample gives a smaller value, v rms — 356±37km/s. 
In our paper we compute the three-dimensional rms pecu- 
liar velocity on scale 8/i _1 Mpc (typical virial radius of a 
galaxy cluster) using Eq. || for the three cosmological mod- 
els. This is because large scale velocities are mostly sensi- 
tive to long wavelength density fluctuations. This smooth- 
ing allows us to get rid of the nonlinear effects on small 
scales but it also tends to underestimate the peculiar ve- 
locities of the smallest objects that we are interested in. 
Nevertheless, with regard to the rather important disper- 
sion in the observational values (320 < v rms < 780 km/s), 
we use the predicted theoretical values, which range be- 
tween 400 and 500 km/s, and are hence in general agree- 
ment with the observational data. 

3.3. Simulations 

For each cosmological model, we generate a simulated map 
of the moving lens effect in order to analyse the contribu- 
tion to the signal in terms of temperature fluctuations. 
The simulations are essentially based on the studies of 
Aghanim et al. (1997). In the following, we describe briefly 
the main hypothesis that we make in simulating the maps 
of the temperature fluctuations induced by the moving 
lens effect associated with small groups and clusters of 
galaxies (10 13 and 10 15 Mq). The predicted number of 
massive objects is derived from a distribution of sources 
using the Press-Schechter formalism normalised (Viana 
& Liddle 1996) using the X-ray temperature distribution 
function derived from Henry & Arnaud (1991) data. This 
normalisation has also been used by Mathiesen & Evrard 

(1997) for the ROSAT Brightest Clusters Sample com- 
piled by Ebeling et al. (1997). The position and direc- 
tion of motion of each object are random. Their peculiar 
velocities are also random within an assumed Gaussian 
distribution. Here again, the correlations were neglected 
because the effect is maximum very close to the central 
part of the structure (about 100 kpc) whereas the correla- 
tion length is between 5 and 20 Mpc (Bahcall 1988). The 
final maps account for the cumulative effect of the mov- 
ing lenses with redshifts lower than z = 1.5. We refer the 
reader to Aghanim et al. (1997) for a detailed description 
of the simulation. 

In this paper, some changes and improvements have 
been made to our previous study (Aghanim et al. 1997). 
In this paper, the predicted source counts (Eq. ||, Sect. 
3.1) are in agreement with more recent data. They are 
also adapted to the various cosmological models that we 
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have assumed. The standard deviation of the peculiar ve- 
locity distribution is computed using equation [| and is 
in reasonable agreement with the data. The advantage of 
using this equation is that the variations with time and 
cosmology are directly handled in the expression. As we 
pointed out in section 2, the secondary effects we study 
here are associated with the whole mass of the structure, 
not only the gas mass. Therefore, the gas part of struc- 
tures are modelled using the /3-profile (as in the previous 
case) to simulate the SZ effect. Whereas the density pro- 
file (Eq. H) is used to simulate the potential well of the 
moving lens effect. We note that the results of the N-body 
simulations of Navarrro, Frenk & White (1996) are con- 
sistent with the assumption of an intra-cluster isothermal 
gas in hydrostatic equilibrium with a NFW halo. 

4. Results of the data analysis 

We analyse the simulated maps of secondary fluctuations 
due to the moving lens effect, for the three cosmological 
models described in Sect. 3, and we quantify their con- 
tributions. We also make attempts at detecting and ex- 
tracting the secondary fluctuations from the entire signal 
(primary CMB, SZ kinetic effect and moving lenses). 



4--1- Statistical analysis 

We show the histogram of the secondary fluctuations 
for the moving lens effect (randomly generated) in the 
three cosmogonies (Fig. ||). In all cases, the amplitude of 
anisotropies ranges roughly between 8T/T ~ — 1.5 10 — 5 
and ST/T ~ 1.5 10 -5 . The rms value of the anisotropies 
varies a little with the cosmological model (43^)^ 5 M — 
5.2 10- 7 , (f ~ 3.4 10- 7 and finally (f ~ 

3.5 10 7 . Our results are in general agreement with those 
of Tuluie, Laguna & Anninos (1996). In all the cosmolog- 
ical models, the rms value of the anisotropies is about a 
factor 10 lower than the rms amplitude of the fluctuations 
due to the SZ kinetic effect associated with the same struc- 
tures, which is about 5. 10 -6 ; and is about 30 times lower 
than the (43-) of the primary fluctuations in a stan- 

V 1 / rms ^ 

dard CDM model. The distribution of the temperature 
fluctuations induced by moving lenses exhibits a highly 
non Gaussian signature (Fig. ||). The fourth moment of 
the distribution, called the kurtosis, measures the peaked- 
ness or flatness of the distribution relative to the normal 
one. We find that the kurtosis for the standard CDM, 
OCDM and ACDM models are positive and respectively 
equal to about 51, 97 and 41. The distributions are thus 
peaked (leptokurtics). 

In the context of our statistical analysis of the sec- 
ondary anisotropies, we also compute the fitted angu- 
lar power spectra (Fig. ||) of the three main sources of 
anisotropies: primary CMB fluctuations (in the standard 
CDM model) and both the predicted power spectra of the 
fluctuations due to the moving lenses (thin lines) and the 
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Fig. 2. Histograms showing the distributions of the secondary 
fluctuations in the simulated maps. The solid, dashed and dot- 
ted lines are for respectively the standard, Open and Lambda 
CDM model. 
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Fig. 3. Power spectra of the primary fluctuations obtained us- 
ing the CMBFAST code compared to the fitted power spectra of 
the secondary fluctuations due to the Sunyaev-Zel'dovich ki- 
netic effect (thick lines) and to the moving lens effect (thin 
lines). The power spectra for the standard CDM model (solid 
line), open CDM model (dashed line) and lambda CDM model 
(dotted line) are shown. 



SZ kinetic effect (thick lines). In figure ||[ the solid lines are 
for the standard CDM model, dashed and dotted lines are 
respectively for the open and non zero cosmological con- 
stant models. We fit the power spectra of the secondary 
anisotropies due to moving lenses with the general expres- 



1(1 + l)Cj = ai s - bis exp(— cut), 



(7) 



in which the fitting parameters for every cosmological 
model are given in table [l]. The SZ kinetic anisotropies 
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Table 1. Fitting parameters for the power spectrum of the 
fluctuations induced by moving lenses as a function of the cos- 
mological model. 





ai s 


bis 




SCDM 


4.910"" 


5.510"" 


3.3 1(T 3 


OCDM 


1.810"" 


7.610"" 


4.910"^ 


ACDM 


2.1 10"" 


2.2 10"" 


3.4 10" 3 



Table 2. Fitting parameters for the power spectrum of the 
fluctuations induced by the Sunyaev-Zel'dovich kinetic effect 
as a function of the cosmological model. 





asz 


bsz 


SCDM 


3.410" 15 


4.310"" 


OCDM 


2.310" 15 


6.310"" 


ACDM 


2.610"" 


2.510"" 



are fitted with the following expression: 

l(l + l)Ci = a sz I + b sz I 2 , (8) 

with the fitting parameters for the cosmological models 
gathered in table. ||. 

The power spectra of the SZ kinetic effect exhibit 
the characteristic I 2 dependence on small angular scales 
for the point-like source dominated signal. All the power 
spectra have rather similar amplitudes, at large scales, 
in particular up to I ~ 200 where we notice an excess of 
power at small angular scales in the OCDM model. This is 
because low f2o models produce higher counts than f^o = 1 
models (Barbosa et al. 1996). 

The moving lens power spectra, for both CDM and 
ACDM models, exhibit a plateau at I > 500 with a de- 
crease at larger angular scales. For the OCDM model, the 
dependence is roughly constant at all scales. We also note 
that the highest and lowest power are obtained, at small 
angular scales, for respectively the standard CDM and 
OCDM models. At large scales, the opposite is true. 
In order to interpret this behaviour, we distinguish be- 
tween what we refer to as the resolved and unresolved 
structures. The spatial extent of the resolved structures 
is much greater than the pixel size (or analogously the 
beam size). Whereas, the unresolved objects have extents 
close to, or smaller than, the pixel size. At the pixel size 
an unresolved structure generates a SZ kinetic anisotropy 
which is averaged to a non-zero value. Whereas the dipo- 
lar anisotropy induced by the moving lens effect is aver- 
aged to zero (except what remains from the side effects). 
A pixel size anisotropy thus does not contribute to the 
signal in the moving lens effect; while it contributes with 
its 8T/T amplitude in the SZ kinetic effect. As a result, 
the distribution of the moving lens anisotropies does not 
reflect the whole population of objects, but only the dis- 
tribution of the resolved ones. In the OCDM model the 



structures are more numerous and form earlier than in a 
standard CDM model. Consequently, the distribution of 
unresolved objects in OCDM thus shows a large excess 
compared with the standard CDM and there are less re- 
solved structures in the OCDM model than in the CDM. 
The excess of power in the the moving lens fluctuations 
spectrum (Fig. [| solid line) reflects the dependence of the 
size distribution upon the cosmological model. 
At a given large scale and for the SZ kinetic effect, there 
is more power on large scales in a standard CDM model 
compared with the OCDM. This is because the contribu- 
tion to the power comes from low redshift resolved struc- 
tures, which are less numerous in an OCDM model. Con- 
sequently, in the case of the fluctuations induced by the 
moving lens effect at large scale, the power in the OCDM 
model is greater than in the standard CDM. In addition, 
at a given large scale the power of the moving lens effect 
accounts for the cumulative contribution from the massive 
objects, with high amplitude, and from the less massive 
ones, with lower amplitudes. 

A comparison between the CMB and the moving lens 
power spectra obviously shows that primary CMB fluctu- 
ations dominate at all scales larger than the cut-off scale, 
whatever the cosmological model (Fig. ||). Furthermore 
in the OCDM and ACDM models the cut-off is shifted 
towards smaller angular scales making the CMB the dom- 
inant contribution over a larger range of scales. The most 
favourable configuration to study and analyse the fluctua- 
tions is therefore the CDM model since it gives the largest 
cut-off scale compared to the other cosmological models 
and since it gives the highest prediction for the power of 
the moving lens effect. The level of spurious additional 
signal associated with the moving lens effect is negligible 
compared to both the primary and SZ kinetic fluctuations. 
Below the scale of the cut-off in the CMB power spectrum, 
the I 2 dependence of the SZ fluctuations is dominant over 
the moving lens effect. Moreover, contrary to the thermal 
effect, the SZ kinetic, moving lens and primary fluctua- 
tions have black body spectra. This makes the spectral 
confusion between them a crucial problem. At small an- 
gular scales, the SZ kinetic effect represents the principal 
source of confusion. 

Nevertheless, the contribution of the SZ kinetic effect 
is very dependent on the predicted number of structures 
that show a gas component. In other words, some objects 
like small groups of galaxies may not have a gas compo- 
nent, and therefore no SZ thermal or kinetic anisotropy is 
generated, but they still exhibit the anisotropy associated 
with their motion across the sky. We attempt to study a 
rather wide range of models. We therefore use two pre- 
scriptions to discriminate between "gaseous" objects and 
"non gaseous" ones. These prescriptions correspond to ar- 
bitrary limits on the masses of the structures. Namely: in 
the first model, we assume that all the dark matter ha- 
los with masses greater than 10 13 Mq have a gas fraction 
of 20% and exhibit SZ thermal and kinetic anisotropies; 
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Fig. 4. Randomly generated power spectra of the primary fluc- 
tuations in the standard CDM model (solid line) compared 
to the secondary fluctuations due to the moving lens effect 
(dashed- dotted line) and the Sunyaev-Zel'dovich kinetic effect 
with a cut-off at 10 13 Mq (dashed line) and at 10 14 Mq (dotted 
line). 

while in the second model, it is only the structures with 
masses M > 10 14 Mq which produce SZ anisotropies. We 
ran the simulations with both assumptions in the stan- 
dard CDM model and computed the corresponding power 
spectra (Fig. |]). The power spectrum associated with the 
SZ kinetic effect shows, as expected, that the cut-off in 
masses induces a decrease in the power of the SZ kinetic 
effect on all scales, and in particular on very small scales 
with a cut-off at I ~ 4000. The power spectrum of the 
SZ kinetic anisotropies can be fitted with the following 
expression: 

1(1 + l)Q = -3.3 10~ 13 + 1.6 10 _14 Z exp(6.2 10~ 4 0- (9) 

Despite this cut-off in mass and the decrease in power, the 
SZ kinetic effect remains much larger than the moving lens 
effect. Therefore at small angular scales, the SZ kinetic 
point like sources are still the major source of confusion. 
In order to get rid of this pollution in an effective way, one 
would need a very sharp but unrealistic cut-off in mass. 

4-2. Detection and extraction 

We analyse the simulated maps in order to estimate the 
amplitudes of the anisotropies associated with each indi- 
vidual moving structure. In such an analysis both primary 
CMB and SZ kinetic fluctuations represent spurious sig- 
nals with regards to the moving lens. Figure || shows that 
these signals contribute at different scales and at differ- 
ent levels. The primary CMB contribution vanishes on 
scales lower than the cut-off whereas the SZ kinetic con- 
tribution shows up at all scales and its power increases 
as Z 2 on small scales. This indicates clearly that the most 
important problem with the analysis of the maps (extrac- 




100 200 300 400 500 

Fig. 5. Plot showing the contours superimposed over a sim- 
ulated map of the fluctuations induced by moving lenses 
in the standard CDM model (pixel size— 1.5'). The 
contour levels and grey scales shown in the plot are: 
1.5 10~ 5 , ±1. 10~ 5 , ±5. KT 6 , ±1. 10" 6 

tion and detection of the moving lens anisotropy) is the 
confusion due to the point-like sources. This problem is 
made worse by spectral confusion. A compromise must be 
found between investigating scales smaller than the CMB 
cut-off, which maximises the pollution due to SZ kinetic 
effect, and exploring larger scales where the SZ contri- 
bution is low (but still 10 times larger than the moving 
lenses) . The main problem here is that on these scales the 
primary fluctuations are 100 times larger than the moving 
lenses which makes their detection hopeless. 

Nevertheless, the signal has two characteristics that 
make the attempts at detection worthy at small scales. 
The first advantage is that the anisotropy induced by a 
moving lens exhibits a particular spatial signature which 
is seen as the dipole-like patterns shown in figure ||. The 
second, and main advantage is that we know the position 
of the center of the structures thanks to the SZ thermal 
effect. 

In fact, the objects giving rise to a dipole-like anisotropy 
are either small groups or clusters of galaxies with hot 
ionised gas which also exhibit SZ thermal distortions. The 
latter, characterised by the so-called Comptonisation pa- 
rameter y, have a very specific spectral signature. It is 
therefore rather easy to determine the position of the cen- 
ter of a structure assuming that it corresponds to the 
maximum value of the y parameter. In the context of the 
Planck multi-wavelength experiment for CMB observa- 
tions, it was shown (Aghanim et al. 1997) that the loca- 
tion of massive clusters will be well known because of the 
presence of the SZ thermal effect. 
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We based our detection strategy for the moving lens effect 
on these two properties (spatial signature and known lo- 
cation). We also assumed that the SZ thermal effect was 
perfectly separated from the other contributions thanks 
to the spectral signature. The problem is therefore eased 
since it lies in the separation of moving lens, SZ kinetic 
and primary CMB anisotropics at known positions. Never- 
theless the clusters and their gravitational potential wells 
are likely to be non-spherical, making the separation diffi- 
cult. In the following, we will show that even in the simple 
spherical model we adopt the separation remains very dif- 
ficult because of the spectral confusion of the moving lens, 
SZ kinetic and primary CMB fluctuations. Separation is 
even more difficult because of the numerous point-like SZ 
kinetic sources corresponding to weak clusters and small 
groups of galaxies for which we do not observe the SZ 
thermal effect. 



4.2.1. Method 

In order to clean the maps from the noise (SZ kinetic and 
CMB fluctuations), we filter them using a wavelet trans- 
form. Wavelet transforms have received significant atten- 
tion recently due to their suitability for a number of im- 
portant signal and image processing tasks. The principle 
behind the wavelet transform, as described by Grossmann 
& Morlet (1984), Daubechies (1988) and Mallat (1989) 
is to hierarchically decompose an input image into a se- 
ries of successively lower resolution reference images and 
associated detail images. At each level, the reference im- 
age and detail image contain the information needed to 
reconstruct the reference image at the next higher resolu- 
tion level. So, what makes the wavelet transform interest- 
ing in image processing is that, unlike Fourier transform, 
wavelets are quite localised in space. Simultaneously, like 
the Fourier transform, wavelets are also quite localised 
in frequency, or more precisely, on characteristic scales. 
Therefore, the multi-scale approach provides an elegant 
and powerful framework for our image analysis because 
the features of interest in an image (dipole pattern) are 
generally present at different characteristic scales. Fur- 
thermore, the wavelet transform performs contemporane- 
ously a hierarchical analysis in both the space and fre- 
quency domains. 

The maps are decomposed in terms of a wavelet basis 
that has the best impulse response and lowest shift vari- 
ance among a set of wavelets that were tested for image 
compression (Villasenor et al. 1995). These two character- 
istics are important if we want to identify the locations 
and the amplitudes of the moving lenses. Since the mov- 
ing lenses induce very small scale anisotropics compared 
to the CMB, we filter the largest scales in order to sepa- 
rate these two contributions. We note that this also allows 
us to separate the contributions due to the large scale SZ 
kinetic sources. In the following we describe our analysis 



method, first applied to an unrealistic study case and then 
to a realistic case. 
Study case 

We filter the large scales of a map of CMB+moving 
lens fluctuations (no SZ kinetic contribution) in order to 
test the robustness and efficiency of the wavelet trans- 
form filtering. In this case, the noise due to the CMB is 
efficiently cleaned. In fact, Figure || lower left panel shows 
a residual signal (symbolised by the dots) associated with 
the moving lens fluctuations, which are simulated in the 
upper right panel of the same figure. We have confirmed 
that the positions of the residual signal agree with the 
positions of the input structures. Moreover, we were able 
to successfully extract the secondary fluctuations due to 
the moving lenses, as well as estimate their average peak 
to peak ST/T values. Figure || shows the average peak to 
peak amplitudes of the input simulated fluctuations (solid 
line) and the extracted values (dashed line). The main 
features are well-recovered, although the amplitudes suf- 
fer from the smoothing of the filtering procedure. In this 
study case, with no SZ kinetic contribution, we find a cor- 
relation coefficient between input and recovered values of 
about 0.95. 
Realistic case 

When this method is applied to filter a map contain- 
ing all contributions (CMB + SZ kinetic + moving lenses), 
we are no longer able to identify or locate the moving lens 
fluctuations, as shown in fig.§ lower right panel. Here, the 
CMB which dominates at large scales is cleaned, whereas 
the SZ kinetic effect, which is mainly a point-like domi- 
nated signal, at least one order of magnitude larger than 
the power of the moving lenses, is not cleaned and remains 
in the filtered signal. We have filtered at several angular- 
scales without any positive result. On large scales the ex- 
tended dipole patterns are polluted by the CMB, as men- 
tioned above, and on small scales the SZ kinetic fluctua- 
tions are of the same scale as the moving lens anisotropies. 
We also tried the convolution of the total map (CMB + SZ 
kinetic + moving lenses) with the dipole pattern function 
but we were still unable to recover the moving lens fluctu- 
ations. In fact, the combination of two SZ kinetic sources, 
one coming forward and the other going backward, mimics 
a dipole-like pattern. In order to distinguish between an 
intrinsic dipole due to a lens and a coincidence, one needs 
to know a priori the direction of the motion which is of 
course not possible. During our analysis, we investigated 
two cases for the cut off in mass as describe in Sect 4.1. 
For the simulations with cut-off mass 10 14 Mq the result- 
ing background due to point-like SZ fluctuations is lower 
than the cut-off at 10 13 Mq case; but we were still unable 
to recover the moving lens fluctuations. 

In our attempt at taking advantage of the spatial sig- 
nature of the moving lens fluctuations, we have located the 
coefficients in the wavelet decomposition that are princi- 
pally associated with the moving lenses and selected them 
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Fig. 6. The two upper panels show simulated maps: on the left, total map of the fluctuations (CMB+SZ+lens) linear scale 
between 1.4 10~ 4 (white) and —1.3 10~ 4 (black). On the right is shown the map of the moving lenses fluctuations (linear scale 
between 1.5 10 -5 (white) and — 1.2 10 -5 (black)). The two lower panels are the result of the wavelet filtering process. The left 
panel is for a CMB+lens configuration, in which we notice that the secondary anisotropics are rather well extracted. The right 
panel represents the case of all contributors and shows that the moving lenses are completely dominated by the SZ kinetic noise. 



from all the wavelet coefficients. Our study case procedure 
is the following. We make the wavelet transform for the 
moving lens fluctuations and, separately, we also make the 
transform for the remaining signals (CMB + SZ kinetic). 
We locate the wavelet coefficients for the moving lenses 
whose absolute values are higher than the absolute val- 
ues of CMB+SZ kinetic coefficients. Then, we select, in 
the transform of the total fluctuation map (CMB + SZ 
+ lenses), the coefficients corresponding to the previously 
located ones. Finally we perform the inverse transform on 
the map (CMB + SZ + lenses) according to the selected 
coefficients. When we compare the average peak to peak 
amplitudes of the recovered (Fig. [?]: dashed line and Fig. 
||: dotted-dashed line) and input (Fig. [j] and Fig. [|, solid 
line) lens fluctuations, we find a very good correlation be- 
tween the amplitudes of the original and the reconstructed 
moving lens fluctuations. The correlation factor is of the 
order of 0.7 for CMB + SZ + lenses with a cut-off mass at 



10 13 M Q and higher than 0.9 with the cut-off at 10 14 Mq. 
This difference between the correlation factors is an effect 
of the cut-off in masses. In fact, for the 10 14 Mq cut-off, 
the filtered maps are cleaner than for the 10 13 Mq cut- 
off. Therefore, in the latter case some of the lenses have 
very little or no signature in the wavelet decomposition, 
hence they are not recovered and the correlation factor 
decreases. 

The results of our study case confirm that the moving 
lens fluctuations have a significant spatial signature in the 
total signal although their amplitudes are very low com- 
pared with the CMB and SZ fluctuations. However, it is 
worth noting that such a "good" result is obtained only be- 
cause we use sorted coefficients from two separated maps, 
one containing the lens signal and the other containing the 
polluting signals. In a real case, there is no way to sepa- 
rate the contributions because of the spectral confusion 
and therefore there is no a priori knowledge of the "right" 
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Fig. 7. Average peak to peak amplitude of the secondary 
anisotropies due to the moving lenses (cut-off mass 10 14 Mq) 
for lenses with a decreasing y parameter. The solid line rep- 
resents the amplitudes in the original simulated lenses map. 
The dashed line represents the extracted amplitudes after sort- 
ing the wavelet coefficients and filtering all contributions 
(CMB+SZ+lenses). The correlation factor is equal to 0.9. 
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Fig. 8. Average peak to peak amplitude of the secondary 
anisotropics due to the moving lenses ( cut-off mass 10 13 Mq ) 
for lenses with decreasing y parameter. Solid line represents the 
amplitudes in the original simulated lens map. The dashed line 
represents the extracted amplitudes without sorting the wavelet 
coefficients and after filtering the CBM+lenses contributions 
(no SZ kinetic). The dotted-dashed line represents the extracted 
amplitudes after sorting the wavelet coefficients and filtering all 
contributions (CBM+lenses+SZ kinetic). 

coefficients in the wavelet decomposition. In our analysis, 
we tried several sorting criteria for the coefficients but we 
could not find a robust and trustworthy criterion to repro- 
ducibly discriminate between the wavelet coefficients be- 
longing to the moving lens fluctuations and the coefficients 
belonging to the noise (SZ kinetic and CMB fluctuations). 



During the analysis, we could not overcome the physical 
limitation corresponding to the presence of sources of SZ 
kinetic anisotropies at the same scale and with amplitudes 
at least 10 times higher than the signal (moving lens fluc- 
tuations). 

5. Conclusions 

In our work, we investigate the secondary fluctuations in- 
duced by moving lenses with masses ranging from those 
of groups of galaxies to those of clusters of galaxies in a 
simple way, based on predicted structure counts and sim- 
ulated maps. This method allows us to explore a rather 
wide range of scales (> 10 arcseconds) in various cosmo- 
logical models. The analysis, in terms of angular power 
spectra, show the scales for which the primary fluctua- 
tions are dominant (Fig. ^|). In the standard and lambda 
CDM models, the primary anisotropies are dominant re- 
spectively for scales I < 4000 and I < 4500 whereas in 
the Open CDM model they are dominant for I < 6000. 
In practice, it is thus impossible to detect the secondary 
anisotropies due to moving lenses in the open model. The 
standard CDM model shows the smallest cut-off scale with 
an intermediate SZ kinetic pollution, compared to the 
other two models. It is therefore the "best case" frame- 
work for making an analysis and predicting the detection 
of fluctuations and the contributions that they induce. 
One must keep in mind that the results quoted in this 
particular case represent the "best" results we get from 
the analysis. 

The results of our analysis are obtained under the as- 
sumption of a universe that never re- ionises, which is of 
course not the case. The re-ionisation, if it is homoge- 
neous, is supposed to somewhat ease the task of extrac- 
tion of the pattern. In fact, its main effect is to damp the 
angular power spectrum of the primary anisotropies on 
small scales, shifting the cut-off towards larger scales. In 
this case, the effect of moving lenses dominates over the 
CMB fluctuations, and the SZ kinetic is not as high as it 
is on very small scales. However, if the re-ionisation is late 
and inhomogeneous, it generates additional SZ kinetic- 
type secondary fluctuations (Aghanim et al. 1996) without 
damping the power spectrum by more than a few percent. 
Here, the re-ionisation might worsen the analysis at small 
scales. In any case, there could be some other additional 
secondary fluctuations principally due to the Vishniac ef- 
fect, that arise in a re-ionised universe. Our work thus 
gives a "best case" configuration of the problem, with all 
other effects tending to worsen the situation. 

We found that the secondary fluctuations induced by 
the moving gravitational lenses can be as high as 1.5 10~ 5 ; 
with rms contributions of about 5 to 3. 10~ 7 in the three 
cosmological models. Even if the moving lens fluctuations 
have a particular dipolar pattern and even if they are "per- 
fectly" located through their SZ thermal effect, the detec- 
tion of the moving lens effect and its separation from the 



12 



N. Aghanim et al.: Moving gravitational lenses: imprints on the CMB 



SZ kinetic and primary fluctuations are very difficult be- 
cause of the very high level of confusion, on the scales of 
interest, with the point-like SZ kinetic anisotropies and 
because of spectral confusion. 

We nevertheless analysed the simulated maps using an 
adapted wavelet technique in order to extract the moving 
lens fluctuations. We conclude that the contribution of the 
secondary anisotropies due to the moving lenses is thus 
negligible whatever the cosmological model. Therefore it 
will not affect the future CMB measurements except as 
a background contribution. We have highlighted the fact 
that the moving lens fluctuations have a very significant 
spatial signature but we did not succeed in separating this 
contribution from the other signals. 
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